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Abstract. We study a Robin boundary problem for degenerate parabolic equation. We suggest a 
notion of entropy solution and propose a result of existence and uniqueness. Numerical simulations 
illustrate some aspects of solution behaviour. Monodimensional experiments are presented. 
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1. Introduction 

Let fi be an open bounded domain of M. N with a Lipschitz boundary and r\ the unit normal to 
<9fi outward to J7. The purpose of this paper is to discuss existence and uniqueness of entropy solution 
for the following initial boundary value problem 

( u t + dwf(u) - A(f>(u) = in Q=]0,T[xfi, 

(P)< u(0,x) =u Q (x) in fl, 

{ b(u)-(f(u)-V<f>(u)).r) = on £=]0,T[x<3fi. 

Here, uq is taking values on [0, M max ] for some w max > 0. Further, the function / is a Lipschitz 
continuous function. Moreover, we require that 

/(0) = and 6(0) = 0. (HI) 

The diffusion term (j) is a continuous function. We consider that there exist a critical value u c of the 
unkown u such that: <f>{.) is zero on [0, u c ] with < u c < u max and (/>(.) is strictly increasing on 
[u c , Umax]- Then problem (P) degenerates to hyperbolic when u takes values in the region [0, u c ] where 
(j) is flat. 

We suppose that the function & is a continuous non-decreasing function on E. In some situation, b 
may be a maximal monotone graph on K (see [4]). Here, we assume also that b satisfies the following 
hypotheses: 

b — ft o (j) where f3 is a non-decreasing Lipschitz continuous function. (H2) 

frOmax) > \f(u max ).ri\- (H3) 
For more than a few decades, the degenerate parabolic equation in bounded domain was studied by 
many authors mainly in the case of Dirichlct boundary conditions (see e.g. [10], [8]). The zero-flux 
boundary condition is studied in [1] for non-degenerate parabolic case, in [7] for fully degenerate 
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hyperbolic equation and recently in [2] for the parabolic-hyperbolic problem. Remark, that the con- 
dition b(u) — (f(u) — V<j){u)).rj = on E includes in particular Neumann (zero-flux) condition on the 
boundary. 

We propose an adequate entropy formulation for problem (P) which incorporates two boundary 
integrals. In [2], existence and uniqueness for the zero flux boundary condition were proved, under 
the assumption (H3) that reads / (tt max ) = in the zero-flux case b = 0. In constrast to the entropy 
formulation in [2], where the passage to the limit in the only boundary integral is straighforward, for 
our entropy inequality, we need the assumption (H2), which permits to give a sense to the boundary 
integral with the term b(u). Indeed, we can deduce that b(u) has a trace on the boundary as a function 
in Sobolev space 

The proof of existence of our entropy solution for any space dimensions N > 1 employs a 
vanishing viscosity approximation. We pass to limit in the interior of the domain Q, by using the local 
compacity result of Panov [12], for this we suppose some relation between / and <p (see Definition 3.5). 
One can refer to [2] for more details. We pay particular attention to the boundary term (here (H2) is 
needed) . 

For the uniqueness result, we use nonlinear semigroup techniques (see, e.g., [6]) and Kruzhkov 
doubling of variables methods. The main goal is to compare two solutions of (P), and it turns out 
that it is simpler to compare a solution of (P) with a regular solution (in the sense that the total flux 
is continuous up to the boundary) of the stationary problem associated to (P). Then we prove that 
entropy solution of (P) is an integral solution, and we refer to the uniqueness of integral solutions 
granted by the general theory of nonlinear semigroup. Unfortunalely, we are not able to obtain regular 
solution to the stationary problem for any space dimensions, but only in one space dimension. Then, 
we can deduce the uniqueness just now when f2 is a bounded open interval of K. Notice that, for the 
same argument as for the zero-flux boundary condition [2], the problem of uniqueness is still open in 
multiple space dimensions. 

The paper is organized as follows. In the next section, we give our definition of entropy solution 
and state some remarks uselful for the well-possedness. In section 3, we prove existence result of entropy 
solution. In the section 4, we prove uniqueness in the case of one space dimension. The latter part is 
devoted to the numerical investigation of problem (P) . We adapt the approach of finite volumes in the 
spirit of Vovelle ([11]) to illustrate and interprete some observations in the case where the assumptions 
(H2) and (H3) are absent. Thereby, we justify the importance these assumptions in this paper. 



2. Notion of Entropy Solution 

Consider the following notion. 



Definition 2.1. A measurable function u taking values on [0, u max ] is called entropy solution of problem 
(P) if <j){u) € L 2 (0,T;H 1 (fl)), b(u) e L 2 (0, T; H 1 ^)) and the following conditions hold: 
Vfc e [0, u max ], V£ e Cg°([0 ) T[xR JV ), with £ > 0: 



Jo Jn I 





+ / |«o- k\£(0,x)dx + / / |/(fc).77(a:) -b{k)\^(t,x)d'H N ~ 1 dt 



(1) 



Here T-L represents the (N — 1)— dimensional Hausdorff measure on dft. 



Remark 1. 1. The entropy solution in the sense of Definition 2.1 is in particular a weak solution. 
Indeed, first take in inequality (1), k = and use (HI). Next, take k = w max and use (H3). 
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2. Let us stress that, in particular, the boundary condition (f(u) — V(/»(u)).r/ = b{u) is verified 
literally in the weak sense as in the case of zero flux boundary condition (see [2]). This contrasts 
with the properties of the Dirichlet problem (see [5]); we expect that the boundary condition 
should be relaxed if assumption (H3) is dropped (see [4, 3]) and numerical tests of section 5. 

3. The integral in the boundary term is well defined due to the hypothesis (H2). We can use the 
fact that the trace of b(u)(t, .) £ on dfl is well defined in L 2 (dVt) for t £ (0,T) a.e. 

According to the idea of J. Carrillo (cf [8]), we give an additional property of entropy solutions, 
useful for the uniqueness techniques. 

Proposition 1. Let £ £ C°°([0, T[xU. N ); then for all k £ [</> c , M max ]; for all D £ and for all entropy 
solution u of (P), we have: 

[ [ {\u-k\£ t + sign(u-k)(f(u)- f(k)-V<p{u)+D).V£}dxdt 
Jo Jn 

+ f \u -k\£(0,x)dx- [ f \b(u)-b{k)\^(t,x)dn N - 1 dt 

Jn JO JdQ 

+ f [ \b(k) - (f(k) - D). v (x)\ttt,x)dH N - 1 dt 
Jo Jan 

> Urn - / / V(t>(u).(V(j>(u)-D)£(t,x). (2) 

In general, uniqueness for evolution equation of kind (P) appear very difficult mainly for the 
initial boundary values problems. In this context, the use of nonlinear semigroup techniques offers 
many advantages. Let us present briefly another notion of solution coming from the theory of nonlinear 
semigroups (see, e.g., [6]). 

Definition 2.2. Let A be an m-accretive operator (see, e.g., [6]). Suppose that h S L 1 (Q), uo & L 1 (fi). 
A measurable function v £ C([0, T]; L 1 (fJ; [0,-Umax])) 1 is an integral solution of the abstract problem 

v t + A(v) 3 h(t), v(t = Q)=u , (3) 

if v(0, .) = u Q (.) and for all (u, z) £ A 

-|K*) - u|| £ i(n) < / sign {v(t) - u)(h(t) - z) + f \h(t) - z\ in V'(0,T). 
at Jn J{v=u} 

We will see that entropy and integral solution coincide in the case = (a, b) an interval of R. 



3. Existence of Entropy solution 

The main result of this part is the following: 

Theorem 3.1. Let N > 1. Assume that (HI), (H2) and (H3) holds. Suppose that (f,(f>) is non- 
degenerate (in the sense of Definition 3.5 below). Then, there exists an entropy solution u for the 
problem (P). 

To show the existence of entropy solutions, we approximate <f)(u) by (j) e {u e ) = <^>(it e ) + eLd(u e ) for 
each e > and set 6 e (M e ) = ft o <j> e (u e ). We obtain the following regularized strictly parabolic problem 
(P e ) with unknown u e 

( u% +div/(u e ) - A0 £ (u £ ) =0 in Q =]0,T[xQ, 

{P e )\ u e (0,x) =u8(a;) in 0, 

[ b t {u e ) - (f(u e ) - V</> £ (w £ )).r/ =0 on £=]0,T[xdQ, 

where (uq) c is a sequence of smooth functions that converges to uq a.e and respects the mini- 
mum/maximum values of uq. 

-'^Hcrc, wc will write L 1 (f2; [0, u m ax]) for the set of all mesurable functions from (1 to [0,u ma x]. 
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Definition 3.2. Let uo be a measurable [0, u max ] -valued function. A function u e E L 2 (0,T; H 1 ^)) 
taking values on [0,w max ] is called weak solution of problem (P e ) if : \/6 E L 2 (0, T; n L°°(Q) 
such that 9 t £ L 2 (Q) and 9(T, .) = 0, one has 

i-T f 

{u t 9t + (f(u e )-W(/) >i (u e )).Ve}dxdt+ I u f Q 8(0,x)dx 



o Jn 

T 

N-l 



b t {u e )6dU N - l dt = 0. (4) 

o Jan 

Theorem 3.3. For uq E [0,u max ], assume (HI), (H2) and (H3) hold. Problem (P e ) admits a weak 
solution u e which is also an entropy solution. In particular, we have < u e < u max . In addition, there 
exists C independent on e such that 

||x/iVt/|| L2(Q) <C; (5) 

\\<l>e(Ue)\\L'(0,T;W(Si)) < C \ ( 6 ) 

||6"(«e)IUi(E) < C and [ u t b n t {u t ) < C. (7) 

This result can be proved, e.g., using Galerkin method (cf [2]). 

Lemma 3.4. Assume that the sequence is such that: ||^ r j||i,2(o,T;ff :L (Q)) < C an d — ^ ^ * n 

L 2 (Q). Then 7^ — > 7^ in L 2 (T,), where 7 is the trace operator. 

The proof uses localization to a small neighbourhood of X. 
To prove existence of entropy solution, we assume that the couple (/(■), </>(■)) i s non-degenerate in the 
sense of the following definition: 

Definition 3.5. (Panov [12]). Let cj> be zero on [0,u c ], strictly increasing on [u c ,u max ] and a vector 
/ = (fx, /jy). A couple (/(•), <P(-)) is said to be non-degenerate if, for all £ E M. N \{0}, the functions 

EN 
£j/i(A) are not affine on the non-degenerate sub intervals of [0,it c ]. 
i—l 

Theorem 3.6. (Panov [12]). Assume that (/, <f>) is non degenerate in the sense of Definition 3.5. Suppose 
u e , e > 0, is a sequence such that 

3d > l,VZ,r 61 with I < r 

T l>r (u e ) t + div(/(IV(u e )) - V</>(TV( W £ ))) is pre-compact in W^f{Q). 

Moreover, suppose u e , f(u £ ), </> £ (u £ ) are equi-integrable locally on Q. Then, there exists subsequence 
(u e ) e that converges in L^ oc (Q). 

Proof of Theorem 3.1. (Sketched) The proof of existence of entropy solution uses Theorem 3.6 to 
justify the passage to the limit in Q (for more details, see [2]) and Lemma 3.4 for boundary integral. □ 



4. Uniqueness result of Entropy Solution in one space dimension 

The main result of this section is the following theorem: 

Theorem 4.1. Suppose that Q = (a,b) is a bounded interval ofM., then (P) admits a unique entropy 
solution. 

In order to study uniqueness in the framework of nonlinear semigroup theory, we consider for all 
bounded function g taking values on [0,w max ], the stationary problem (5*) associated to problem (P): 

( u + div(/(u) - V0(u)) =g in n, 
[ ' \ b(u) - (f(u) - V<f)(u)).T] = on dn. 
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The notion of entropy solution of (S) correspond to the time-independent entropy solution of (P) 
with source term g — u. In the case where CI = (a, b) is a bounded interval of R, we have an important 
result, which states that, the total flux is regular at the points a and b. This kind of regularity seem 
hard to obtain in multiple space dimensions for (£), and even in dimension N = 1 for (P). 

Proposition 2. For all measurable function g taking values in [0, M ma x] the problem (S) admits a 
solution u such that (f(u) — 4>(u) y ) is continuous up the boundary, i.e.. (f{u) — 4>{u) y ) G C([a,b]). 
Moreover, b(u) — (f(u) — (j)(u) y ).ri(y) is zero at y = a and y = b. (Here 77(a) = —1 and rj(b) = +1 ). 

From now, let's define the operator Af^ t t, on L 1 associated with regular solutions of (S) by its 
graph: 

(u, z) £ Af,<j>,b ={u such that u is an entropy solution of (S), with g = u + 2}. 

Proposition 3. 1. Af^j, is accretive in i 1 (S7) . 

2. For all A sufficiently small, R(I+ XAf ^ j,") contains L (CI; [0, u max ]). 

3. D(A f ^ b ) = L 1 (Cl;[0,u max }). 

For the proof of this proposition, we can refer to [2]. 
According to the general results of [6], it follows existence and uniqueness of integral solution in the 
sense of Definition 2.1: 

Corollary 1. Let CI — (a,b), uo,uq £ L 1 ^) andh,h G L l (Q). Letv,v be integral solutions of (3) (with 
operator Af t j,^ ) associated with the data (uq, h) and (uq, h), respectively. Then for a.e. t £ [0,T). 

\\v(t) - v(t)\\ L i <\\u - u \\ L i + f \\h(r) - h(T)\\ L idt. 

Jo 

Adapted to our case, we have the following result 

Theorem 4.2. Let CI = (a, b). Let v be an entropy solution of (P) and u be an entropy solution of (S). 
Then 

^||u(t) - u||z,i(n) < J sign(v -u){u- g)dx in V'(0, T). (8) 
In particular, v is an integral solution of (3) with h = 0. 

Proof of Theorem 4.2 and Theorem 4.1. We consider v = v(t,x) an entropy solution of (P) and u = 
u(y) an entropy solution of (S). Consider nonnegative function £ = £(t, x, y) having the property that 
£(.,., y) G C°°([0,T) x O) for each y £ CL, £(t,x,.) £ C^(U) for each (t,a;) G [0,T) x fi. Apply the 
doubling of variables [9] in the spirit of [2], we obtain this following inequality 

m\v - u\£ t dydxdt + / / \v - u\£(0, x, y)dxdy 
.. a JnJn 

sign(v - u) (f(v) - (f>(v) x ) - (f(u) + <j>{u) y ) + £ y ) 



JnJn 

T 

\b(u) - (/(«) - <f>(u) y ).ri(x)\t;dydo-dt 

xeenJn 

\b(v) - (f(v) - <t>(v) x ).ri(y)\Sdo-dxdt 

/q JnJyedn 

+ / / sign(v ~ u)(u — g(y))£dydxdt 
Jo JnJn 

> f f [ \b(u)~b{v)\£do-dxdt+ [[ I \b(u)-b{v)\£do-dydt 
Jo Jxcn Jyedn Jo J yen Jxedn 



lim- / // \<f>(v) x -cj){u) y \ 2 Zdydxdt>0. (9) 

lngxn^n{-<j<c/>(v)~4>(ii)<cr} 



<7->o a 
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Next, following the idea of [1], we take the test function £(t,x,y) = 9(t)p n (x,y), where 8 SC§° (0, T), 
9>0, p n (x,y)=S n (A) and A=(l-^^)x-y+^ s} . Then, Pn G 25(0x0) and Pn lnxsn = 0. 
Due to this choice, 

/ / / |&(«)- C/(«)-^)*)-»7(»)|PnW^* = 0- 
JO JxSfi JyGdfl 

By Proposition 2, f>(tt) — (f(u) — <f>(u) y ).r](y) € Co([a, &]). Therefore we have 
|6(it) — (f(u) — (j)(u) y ).r)(x)\ — > when x — > y, i.e, as n — > oo. We conclude that 

lim f f [ \b{u)-(f(u)-<t>(u) y ).Ti(x)\p n 9dydadt = 0. 

n_>oo J0 JxedQ JyGQ 

with the calculation detailed in [2], we deduce that 

Osign(v - u) (f(v) - <j>(v) x ) - (f(u) - <p{u) y ) .{{p n )x + (p n ) y )dydxdt -> 0. 

! L J 

Hence, we get (8) by passing to the limit in (9) with the above choice of £. Thus, the entropy solution 
v of the problem (P) is an integral solution of (3). This proves that v is a unique entropy solution due 
to Corollary 1. □ 



T 

o JQ.JQ. 



5. Role of Hypotheses (H2), (H3) and some numerical illustrations 

The numerical analysis of (P) is not the aim of this paper, although we consider this alternative in a 
future work. We present briefly the importance of the hypotheses (H2), (H3). We apply now the ideas 
developped e.g., in the work ([11]) of Vovelle to construct a monotone scheme which take into account 
the boundary condition. The interval [0,1] is divided into I cells and the numerical approximation at 
t = nSt in the cell number i is : 

-t pi5x 

Vi € {1, ...,/}: u°i = — / u a (x)dx (10) 



< +i = <- ^ (p«,< +1 )- — 6x 1 ' v % ~ u ) mi 

with the boundary conditions taken into account via 

< +1 = «?- £ (pm, «?)- 0K) ;/ K) - 6 K) ) . (12) 

,., <5t / 0« i) 

f w+l . „.n ( l/-.?^ T*f~.n „.n\ , i / ^ ^ 1 — 1/ 



""'far W) ' fW '""" )+ ^ (13) 



Here, P is a numerical flux which we assume monotone, consistent, Lipschitz continuous (see [11]). 
In the sequel, we take uq(x) = 0.7 if x G [|,1] and uo(x) — if x G [0, |[ for Figure 1, Figure 2, 

Figure 3 and Uo(a;) = sin(y) for Figure 4. We take <j>{u) = (u — 0.6) + , Sx = 0.01, St = ^f, and 
consider a numerical solution at time t = 0.1. Initially, we remove the hypothesis (H3), by taking 

2 

f(u) — \ and b(u) — <fr(u), we observe numerically the loss of maximum principle (see Figure 1 ). 
Our entropy formulation requires to choose b(u) in the functional space that permit to define the 
trace of b(u) on the boundary. In the context where assumption (H2) is not taken into account, 
b(u) = u and f(u) — u(l — u)l[ .i]; numerically, we observe a boundary layer (see Figure 2 ) and this 
is confirmed by theorical results of [4]. Now, taking into account assumptions (H3), (H2), with data 
f(u) — u(l — u)liQ t xy, b{u) — <p(u) the numerical observation shows that the boundary condition at 
x = and x = 1 is verified literally and the numerical solution respect the maximum principle (see 
Figure 3 and Figure 4). 
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Figure 1. 



Figure 2. 



Figure 3. Figure 4. 
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